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Abstract. We investigate the cross-correlation between the cosmic infrared and microwave 
backgrounds (GIB & CMB) anisotropies through the integrated Sachs-Wolfe effect. We first 
describe the CIB anisotropies using a linearly biased power spectrum, then derive the theoretical 
angular power spectrum of the CMB-CIB cross-correlation for different instruments and 
frequencies. We discuss the detectability of the ISW signal by performing a signal-to-noise 
(SNR) analysis with our predicted spectra. The significances obtained range from 6a to 7a in 
an ideal case, depending on the frequency ; in realistic cases which account for the presence of 
noise including astrophysical contaminants, the results span the range 2-5a, depending strongly 
on the major contribution to the noise term. 



1. Introduction 

The discovery of the acceleration of the expansion of the Universe, made through supernova 
observations at the end of the last century, led us to hypothesize the existence of an unknown 
"Dark Energy" , contributing 70% to the energy budget of our Universe. In this presentation, we 
focus on a specific probe of the Dark Energy, namely the integrated Sachs- Wolfe (ISW) effect : 
This effect of gravitational origin induces secondary anisotropies in the CMB and is due to large- 
scale structures. The gravitional potentials of the latter are slowly decaying in a A-dominated 
Universe, and therefore give a net difference in energy to the CMB photons that travel across 
them. This effect shows in the power spectrum of the Cosmic Microwave Background (CMB) 
temperature anisotropies at large angular scales but its small amplitude and the cosmic variance 
at those i make its direct detection very challenging, if not impossible, when using only the CMB 
itself. To circumvent this limitation, cosmologists have devised a way to exploit the link between 
this imprint on the CMB and the large-scale structures causing it, by simply cross-correlating 
the CMB with matter density maps (galaxy maps in practice) and then comparing the results 
to a null hypothesis and to what is expected from theory. 

During the last decade or so, a growing interest has risen in this field thanks to the 
development of large galaxy surveys in many wavelengths. However, this method has yet to 
produce a definitive and conclusive detection of the ISW effect from current surveys, with 
significances ranging so far from negligible up to more than 4a, and with sometimes confiicting 
results throughout the literature. We make here a short report of our work presented in Ilic 
et al., 2011 ([5)- The originality of our work is to consider the Cosmic Infrared Background 
(CIB), first discovered by Puget et al. in 1996. This background, visible roughly from 10 to 
1000 /im in wavelength, arises from the accumulated emission of star-forming galaxies spanning 
a large range of redshifts including our current Dark Energy dominated era. The CIB also 
features anisotropies that are underlined by the galaxy density field and thus the matter density 
fiuctuations. It is therefore reasonable to expect that it has a positive correlation with the CMB 
through the aforementioned ISW effect. 



2. Modelling the expected signal 

2.1. CIB anisotropies 

Ever since its discovery, many efforts have been deployed to detect tlie CIB and its anisotropies 
witli increasing precision. Tlie most recent studies on tlie subject make use of tlie halo occupation 
distribution and the Dark Matter halos properties, in order to predict the power spectrum of 
these anisotropies. Such models are particularly useful when describing the small, non-linear 
scales of the CIB. Since we focus here on the ISW effect which only concerns much larger scales, 
we can use a simpler model for the CIB anisotropies power spectrum, similar to the description 
made by Knox et al. in [2]. We first write the CIB temperature anisotropies as the following 
line-of-sight integral : 



with 5j being the emissivity fluctuations of the CIB. Again, similarly to [2J, we assume 
that the CIB anisotropies are direct tracers of the matter density fluctuations up to a factor 
bj{i', z), a frequency- and redshift-dependent matter-emissivity bias : 5j{{r]o — r/)n, z/, z) z) = 
hj{u, z)6{{riQ—T])h, z) where z) is the mean emissivity per comoving unit volume at frequency 
as a function of redshift z, and is derived using the empirical, parametric model of [3]. We 
describe the matter density field S in our analysis by a linear power spectrum valid for the 
scales of interest in our work (i < 100, comprising most of the ISW signal). Lastly we choose 
the linear bias hj{v,z) to be constant in redshift : bj{i',z) = 6iin(z^)- As detailed in Ilic et al. 
([!]), to obtain it at each frequency, we compute the value of bun that gives the best agreement 
between our predicted linear CIB power spectrum and the one obtained from the Planck data 
and published in [4J. We get the following values for these bias factors : bu^ = 1.74 at 857 GHz, 
2.09 at 545 GHz, 2.63 at 353 GHz and 2.51 at 217 GHz ; these values are consistent with the 
results found in the literature. The linear bias we obtain this way increases with the wavelength : 
this is coherent with the fact that as we go deeper into the infrared, the galaxies probed are 
more luminous at higher z. They reside in more massive and rarer halos, and are therefore more 
biased. 

2.2. Correlation with the ISW 

In the CMB anisotropies, the temperature contribution due to the ISW effect is an integral over 
the conformal time of the growth rate of the gravitational potentials : 



where 7]r is some initial time deep in the radiation era, ^ and ^ are the Newtonian gauge 
gravitational potentials, r(ry) is the optical depth, and the dot denotes differentiation with 
respect to rj. Using this equation and Eq. ([1]), and going to Fourier space, we can express the 
CMB-CIB cross-power spectrum at a frequency u : 



where comes from the primordial curvature power spectrum Pn = 27r^A^/A;^ ; j^(-) are the 
spherical Bessel functions, while S, 4> and '0 are the time-dependent parts of (respectively) the 
matter density contrast 6, and the two Newtonian gravitational potentials $ and ^. The two 
coefficients C\ii$ and Qij, give the relations between 5, <I> and ^ for adiabatic initial conditions : 
c<5* = <5/^ = —3/2, c^r<j) = ^'/<I> = —(1 -I- {2/5)Ru), where Ri, = Pv/{pv + P7), with and p^ 
respectively the energy densities in relativistic neutrinos and photons. 




(1) 





3. Signal-to-Noise analysis 

To compute the aforementionned spectra, we adapted for our analysis a code named 
CROSS_CMBFAST [5] : for a given cosmology and emissivity function j{i',z), our code 
calculates the C^^ from Eq. ([3]) and at the same time the predicted power spectrum of the 
CIB fluctuations described by Eq. ([1]). 



3.1. Ideal case 

We now investigate the detection level of the ISW effect using CMB-CIB cross-correlation by 
performing a signal-to-noise ratio analysis in the context of several past (IRAS) and present 
{Herschel SPIRE, Planck) experiments and their associated frequencies. Using our computed 
power spectra, we can write for each given frequency ly the total signal-to-noise ratio of the ISW 
detection as : 

where the total (or cumulative) signal-to-noise is summed over multipoles between 1 = 2 and 
^max ^ 100 where the signal has its major contribution. 

We first consider the ideal situation where the CIB and CMB maps used for cross-correlation 
are noiseless and cover the whole sky ; with these ideal assumptions, we obtain high levels of 
detection for the CIB-CMB correlation which reach ~ 6-7fT, with the largest contribution to 
the SNR coming from multipoles lower than ~ 50. Our results for the Planck instrument are 
shown on Table [1] ; in light of these, the optimal frequency for ISW detection appears to be 
around 353 GHz with a maximum SNR reaching 7a. However in pratice, the CIB extraction at 
this frequency might prove challenging since the CMB becomes dominant here, and increasingly 
so as we go down in frequency. Therefore the possible residuals in the extracted CIB map have 
to be accounted for, and other sources of noise as well, which we do next. 



3.2. More realistic SNR 

We now carry a more realistic study by including several possible sources of contamination : first, 
the strong emission of the Milky Way masks the CIB signal over the galactic plane, reducing 
the usable fraction of the sky (/sky) by at least ~ 25%, and consequently reducing the SNR by a 
factor y^/sky. Furthermore, the rest of the sky is also quite polluted by dust foregrounds, which 
will have to be removed from our CIB maps although some residuals might remain. There may 
even be a significant CMB residual in this map, due to an imperfect separation of components. 
Consequently, we assess the impact of these contaminants in our study, by adding a noise term 
A'^^^ to the CIB power spectrum in our SNR. We break this noise component into several 
independant parts : 

where these four different terms represent the power spectra of the CMB residual, the 
galactic foregrounds residuals, the instrumental noise and finally the noise due to correlation 
between residuals and the CIB (which appears when autocorrelating the final CIB map). The 
instrumental noise is fixed by the characteristics of the considered experiment. We then choose 
three parameters to govern the other three contributions : /sky, the fraction /cmb of the total 
CMB map present as a residual in the CIB map, and ^fore. the amplitude of the foreground 
residuals spectrum relative to the CIB spectrum (at i = 10). The next step would be to explore 
this 3D parameter space at each frequency and compute the SNR at each point. Considering 
the very large number of possible combinations of parameters, it would not be practical to 
display the complete results of this exploration here. Therefore we first choose to fix /sky to 
two values of interest : /sky = 0.75, corresponding to an optimistic case where only the galactic 
plane is discarded ; /sky = 0.15, a pessimistic estimate of the area of the sky where the current 
data allows for an efficient CIB extraction. For our other two parameters we limit ourselves to 
reasonable values, with /cmb £ [0,0.1] and ^forc. £ [0, 10]. We also explore the combination of 
cross spectra at different frequencies, which allows us to increase the total signal-to-noise ratio of 
the ISW detection by combining their constraints, although this will be limited by the possible 
intrinsic correlations between the different CIB frequencies. 



Table 1. Total signal-to- noise ratio of the CIB-CMB cross-correlation for four of the CIB 
frequencies of Planck HFI and three different noise scenarii. 



Frequency (GHz) 


857 


545 353 


217 


Wavelength {f.im) 


350 


550 850 


1380 


Perfect Single SNR 


6.26 


6.83 6.98 


6.95 


Joint SNR 




7.12 




Realistic SNR 1 (/^ky = 0.75, 


5.36 


5.73 5.39 


3.56 


/cMB = 0.01, Aorc. = 0.01) 








Joint SNR 




5.88 




Realistic SNR 2 (/^ky = 0.15, 


2.40 


2.56 2.41 


1.59 


/cMB = 0.01, Afore. = 0.01) 








Joint SNR 




2.63 





After a thorough exploration of these parameters, we show a part of our results in Table [TJ 
and we can draw a few conclusions : First the gain brought by the joint correlation is small, 
which can be attributed to the high correlation (therefore redundancy) between the CIB at 
its different observed frequencies. When foucusing on the 353 GHz results (the best SNR in 
the ideal case), the influence of the CMB is clearly visible, quickly reducing the SNR as its 
residual level increases, which is even more pronounced at 217 GHz, due to the proximity to the 
maximum of the SED of the CMB. The presence of instrumental noise becomes significant at the 
two lowest frequencies (217 and 353 GHz), again reducing their value in the cross-correlation. 
As expected the galactic foreground residuals also decrease the SNR, though their influence is 
roughly the same at all frequencies as they are deflned relatively to the CIB spectrum. Lastly, 
the biggest influence comes from the fraction of the sky through the /sky parameter, as the total 
SNR scales as /sky This makes it a crucial requirement for future applications to have the 
largest possible coverage to minimize this effect. 

4. Conclusions 

We investigated the cross-correlation between the cosmological infrared and microwave 
backgrounds through the ISW effect, and studied its detectability under various observational 
situations. Using an advanced SNR analysis which included the main sources of noise both 
instrumental and astrophysical, and all their possible correlations, we pointed out the most 
promising frequency in the ideal case of noiseless full-sky maps (353 GHz) with an expected 
signiflcance as high as 7a for the cross-correlation signal. The same frequency turned out to be 
less optimal with more realistic assumptions about sky coverage and possible sources of noise 
(here CMB, dust residuals and instrumental noise). In this case, higher frequencies such as 
Planck HFI's 545 and 857 GHz are favored, with an expected signiflcance ranging from 2 to 
more than 5a depending on the frequency, the levels of noise and the fraction of the sky available 
for analysis. We also found that a joint cross-correlation using all available frequencies is of minor 
interest. Nevertheless, our best results for /sky = 0.75 are higher than the signiflcances of all 
current CMB-galaxies cross-correlation, with a > 5, although a less optimistic estimate for the 
sky coverage quickly reduces our signal-to-noise ratios. This stresses once again the requirement 
of good component separation techniques and foreground removals for future applications, in 
order to have the largest fraction of common clean sky /sky possible. We are currently applying 
our method and framework to the actual Planck data, and the results of this study will be 
published with the next public release of Planck papers and data at the beginning of 2013. 
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